Quantum Control of the Hyperfine Spin of a Cs Atom Ensemble 
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We demonstrate quantum control of a large spin-angular momentum associated with the F = 3 
hyperfine ground state of 133 Cs. A combination of time dependent magnetic fields and a static 
tensor light shift is used to implement near-optimal controls and map a fiducial state to a broad 
range of target states, with yields in the range 0.8-0.9. Squeezed states are produced also by an 
adiabatic scheme that is more robust against errors. Universal control facilitates the encoding and 
manipulation of qubits and qudits in atomic ground states, and may lead to improvement of some 
precision measurements. 
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Accurate dynamical control plays a central role when 
quantum mechanics is leveraged to improve the outcome 
of some physical process. Quantum control has been ac- 
complished in many contexts and at various levels of so- 
phistication, and has become a mainstay in areas such as 
nuclear magnetic resonance coherent chemistry Q, 
quantum information processing Q and quantum metrol- 
ogy 01 • O ne extensively studied problem is how to trans- 
fer a physical system from an initial fiducial state (e.g. 
the ground state) to some final state, as is done, for ex- 
ample, in optical control of chemical reactions Q. In such 
cases the figure of merit for control is the yield, or overlap 
between the actual and desired states. As long as errors 
and decoherence are negligible the general topography of 
control landscapes (yield vs control parameters) is well 
understood [B| and techniques are available for efficient 
design of optimal controls [6(. The most ambitious level 
of quantum control requires that the system be control- 
lable in the Lie-algebraic sense Q, a sufficient condition 
for which is that internal dynamics plus interaction with 
external fields can generate any unitary transformation 
within the Hilbert space of interest. Even when full con- 
trol is possible in principle, attention must be paid to 
robustness in the presence of realistic levels of dissipa- 
tion and systematic errors in the control fields. In spin- 
1/2 systems this can be accomplished through open loop 
control [H, i. e. without recourse to real-time feedback 
Q or error correction but little is known about how 
to systematically design robust controls in a larger state 
space. 

In this letter we demonstrate quantum control of the 
spin-angular momentum (nuclear plus electronic) associ- 
ated with the F = 3 hyperfine ground state of individ- 
ual 133 Cs atoms, corresponding to a 2F + 1 = 7 dimen- 
sional Hilbert space. Starting from an easily prepared 
fiducial state we use magnetic fields and AC Stark shifts 
(light shifts) to design and implement near-optimal con- 
trols and produce a range of target states. We evaluate 
our control performance by experimentally reconstruct- 
ing the entire spin density matrix [l3| and computing the 



overlap between the measured and target states. In most 
cases the estimated yield is in the 0.8 — 0.9 range, limited 
by errors in the control fields and to a lesser extent by 
decoherence from light scattering. The measured states 
can be compared also to the predictions of a full model 
that includes the effects of errors and decoherence. Typ- 
ical fidelities between measured and predicted states are 
around 0.9, which is close to the resolution limit of our 
procedure for quantum state estimation. We further use 
this universal approach to generate spin-squeezed states 
and compare against a method based on adiabatic evolu- 
tion [3] ■ The latter is more robust against errors in the 
control fields, but also slower and thus more sensitive to 
light scattering and decoherence. Large spins provide a 
testing ground for the design of accurate and robust con- 
trols in a system where the Hamiltonian is well known 
and where errors and dissipation are well understood and 
can be accurately modeled. From a practical perspective, 
quantum control of hyperfine states has direct relevance 
to proposals for neutral atom quantum compu ting (loj 
wherein qubits (or higher dimensional qudits [111 ]) are 
encoded in the ground-state manifold, and may provide 
a simple route to modest spin squeezing and accompany- 
ing gains in precision atomic magnetometry [T3 |. 

Universal control of a spin F requires that the Hamil- 
tonian dynamics be capable of generating an arbitrary 
unitary transformation in SU(2F + 1). A linear Zee- 
man interaction between the atomic magnetic moment 
and a moderate strength magnetic field generates only 
Larmor precession and the geometric rotations that are 
a representation of 577(2). More general control requires 
a Hamiltonian that is non-linear in at least one compo- 
nent of F. In our experiment this is provided by an off- 
resonance light field that couples to the atomic ground 
state through the tensor ac-polarizability and leads to 
a spin-dependent light shift with an irreducible rank-2 
component [15]]. The combination of a time-dependent 
magnetic field and a constant x-polarized light field re- 
sults in a control Hamiltonian [16fl, 
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FIG. 1: Quantum control of a large atomic spin, (a) Schematic of the experiment, (b) Example of a control waveform <j>(t). 
(1-4) Wigner functions at four stages during the control sequence. Both sides of the sphere are shown. The final result is close 
to the target state \xt) = (\m z = 2} + jm z = —2})/y/2. (c) Density matrix (absolute values) and Winger function for \xt) 
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where we have expressed the strength of the nonlin- 
earity in terms of the photon scattering rate 7 S and 
where the dimcnsionless parameter (3 is a measure of the 
timescales for coherent versus incoherent evolution. Its 
value depends on the atomic structure and the frequency 
of the driving field and for Cs takes on a maximum value 
(3 = 8.2 when tuned between the hyperfine transitions 
of the D\ line at 894nm. This is enough to allow con- 
siderable coherent manipulation. It follows directly from 
the theory of Lie groups that a Hamiltonian of this form 
allows full control of a spin of any magnitude 17[ • Specif- 
ically, one can show that the algebra generated by com- 
mutators and linear combinations of F Xl F y , F x spans the 
entire (2F+ 1) 2 — 1 dimensional operator space necessary 
to represent SU(2F+ 1). Thus, a time-varying magnetic 
field in the x — y plane suffices to make Hq (~t) universal. 

A schematic of our setup for spin quantum control is 
shown in Fig. QJa). We begin with a sample of a few 
million Cs atoms, captured and laser cooled to ~ 2/z-fT in 
a magneto-optical trap and optical molasses. Once the 
atoms are released from the optical molasses their spin 
state is initialized by optical pumping into a state of max- 
imum projection along the y-axis, |V>o) = \F = 3,m y = 
3). We drive the spins by applying a time-dependent 
magnetic field from a set of low-inductance coils driven by 
arbitrary waveform generators, and by applying a static 
light shift from an optical probe beam. Using an all-glass 
vacuum cell, avoiding nearby conductive or magnetizable 
materials, and synchronizing our ~ 0.5ms duration ex- 
periment to a fixed point during the AC line cycle allows 
us to null the background magnetic field to a few tens 



of /xGauss without the use of shielding or active com- 
pensation. The applied magnetic field can be controlled 
with an accuracy better than one percent in a bandwidth 
of more than 100kHz. Immediately following a period of 
quantum control we estimate the resulting quantum state 
as described in [131 ] . In this procedure the control mag- 
netic and optical fields are applied to drive the spins for 
an additional 1.5 ms, while continually and weakly mea- 
suring a spin observable through its coupling to the probe 
polarization. To reduce the effect of noise, the measure- 
ment signal is averaged over 16 repetitions of the exper- 
iment and the full density matrix determined from the 
measurement record and the known evolution. 

Control Hamiltonians for our experiment are designed 
through a simple procedure that we have found empiri- 
cally to produce very good though not provably optimal 
results. The objective is to start from the state \tpo) and 
to produce a specified target state \xt) by modulating 
the field B(£) for a fixed time r. With readily available 
magnetic fields the timescale for geometric rotations is 
much shorter than for nonlinear evolution driven by the 
light shift, and the latter therefore becomes the time- 
limiting element of most transformations. In our experi- 
ment the maximum available Larmor frequency is 15kHz 
and the nonlinear strength is (3js ~ 2-7T x 500-ffz. Under 
these conditions there is no significant sacrifice in control 
performance when the set of available rotations is some- 
what restricted. We therefore choose the magnetic field 
to have constant magnitude and time-dependent direc- 
tion in the x-y plane. With this simplification the con- 
trol Hamiltonian is completely determined by the time 
dependent angle 4>{t) between B(£) and the s-axis. The 
state \xt) and the transformation Vo) —> \xt) belong to 
a d = 7 dimensional Hilbert space and can be specified 




FIG. 2: Examples of target and measured density matrices 
(absolute values). The target states are (a) {\m z — 2) + \m z — 
—2})/y/2 , (b)|m^ = 2), (c) 'S, y m y \m y ). The experimental 
yield is indicated for each case. 



by a set of 2d — 2 = 12 real numbers, and full control 
therefore requires at least that many free parameters in 
the control Hamiltonian. To ensure sufficient flexibility 
we specify the control waveform <fi(t) by its values {4>i} 
at N = 30 discrete time steps. 

The design of a control waveform now proceeds 
through two search iterations. In the first round we cal- 
culate the state \tpp) produced by a sequence of field 
directions {(pi} by integrating the Schrdinger equation 
for a time r, with suitable filtering of the corresponding 
B(t) to reflect the bandwidth and slew rate limitations 
of our magnetic coils and drivers. A locally optimal con- 
trol waveform is found by starting from a random seed 
and maximizing the yield *W = \(xt\"4>p)\ 2 with a sim- 
ple gradient ascent algorithm. We have found that a 
small set of random seeds almost always generates at 
least one waveform with yield greater than 0.99, which 
is expected from the general structure of control land- 
scapes derived in [Bj]. At this point we switch to a more 
realistic estimate of control performance by modeling the 
evolution with a full master equation that incorporates 
decoherence from light scattering and inhomogeneity of 
the nonlinear strength across the atomic ensemble. This 
allows a second stage of optimization starting from the 
waveform generated in round one and using the more 
complete but computationally intensive model to predict 
the yield, which is now defined in terms of the overlap 
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p T and the predicted density matrix p p . 

An example of an optimized control waveform is shown 
in Fig. |Tfb) , along with Wigner function representations 
of the spin wavepacket jl8| at a few steps during the 
transformation as calculated using the complete master 
equation. Note that the nonlinear evolution initially pro- 
duces a squeezing ellipsoid which later wraps around the 
sphere so that interference effects can be manipulated to 
create the desired state. The end product is very close 
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FIG. 3: Histograms of (a) yields and (b) fidelities of measured 
vs. predicted states, (c) & (d) Yields and fidelities when each 
measured state is geometrically rotated to optimize overlap 
with the predicted state. 



to the target state shown in Fig. [He). According to our 
model this and a wide variety of other control waveforms 
all produce yields near 0.95. Taking into account im- 
perfect optical pumping in our experiment (the initial 
population in \ipo) is ~ 0.96) reduces the expected yields 
to around 0.90. 

We have generated and tested a sample of control 
waveforms designed to produce 21 different pure spin 
states. Fig. [2] shows three examples of target and mea- 
sured density matrices, with yields falling in the range 
0.87-0.97. A more complete statistics of yields for over 
a hundred experimental realizations of control is com- 
piled in the form of a histogram in Fig. [3](a) , show- 
ing a fairly broad distribution centered on respectable 
value of 0.8. It is also informative to compare the mea- 
sured density matrices p M against the density matrices 
p P predicted by our model, as quantified by the fidelity 

Fig. [2Kb) shows a histogram of 



T rp I 1/2 1/2 

T = Tr\/p p ' p M p p ' 



fidelities for our data set. Note that both yield and fi- 
delity can be affected by control errors (the real state 
is different from p P ) as well as state estimation errors 
(the real state is different from p M ), and that there is 
no way to distinguish between these possibilities. Nu- 
merical modeling shows that small background magnetic 
fields or miscalibration of the control fields will lead to 
apparent geometric rotations of the final state, but such 
errors are too small in our experiment to significantly 
affect the outcome. The obvious outliers in the yield 
and fidelity distributions are associated with two specific 
control waveforms, and closer examination shows that 
the estimated states are rotated relative to the predicted 
states. The axis of rotation corresponds the direction of 
the magnetic field at the transition between the control 
and state estimation phases, which suggests a problem 
with the way the corresponding control waveforms were 
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joined together. We can numerically rotate a given p M to 
maximize its fidelity relative to p P and obtain new values 
for yield and fidelity. Carrying out this procedure for all 
data points takes care of the outliers without otherwise 
changing the yield distribution significantly, as shown in 
Fig. [3f c) . This distribution can reasonably be interpreted 
as a measure of our ability to control the spins in a well 
designed experiment. The fidelity distribution (Fig.|3(d)) 
remains peaked at 0.9, which we know from experience 
to reflect the accuracy of our state estimation algorithm. 
Finally we note that random errors in state estimation 
are far more likely to decrease than increase the appar- 
ent yield. A simple error model based on Gaussian ran- 
dom displacements in state space indicate that the yields 
are probably 10% larger on average than indicated by 
Fig. E^d). This puts most yields in the range 0.8-0.9, in 
good agreement with the ~ 0.9 predicted by the model 
used to design the control waveforms in the first place. 

To further explore quantum control in our system we 
have studied the generation of spin squeezing both by 
optimal control as outlined above and by the adiabatic 
scheme described in [3]. The latter begins with an ini- 
tial state, \ipo) — \F = 3,m y — —3), which has equal 
uncertainties for the components AF X and AF Z and is 
often referred to as a spin-coherent state. This state is 
a good approximation to the ground state of the control 
Hamiltonian He (t) when the magnetic field is of the form 
B(i) = B(t)y and B(t) is large. As the field magnitude 
is slowly reduced the state adiabatically evolves so as 
to minimize the squeezing parameter £ = AF x /\(F y )\ of 
relevance for metrology [191 ] . Fig. BJa) shows the pro- 
gression of squeezing and anti-squeezing relative to a 
spin-coherent state with the same |(-Fj,)|. Up to ~ AdB 
of squeezing is seen in the experiment, in good agree- 
ment with the predictions of our model. For the small 
spin magnitude used here the squeezing is quickly limited 
by the decrease in |(Fj,}| as the squeezing ellipse wraps 
around the sphere. Fig. H)^b)-(c) shows Wigner functions 
for the target and actual state for the smallest £ achieved 
in our experiment (~ 80% of the coherent state value). 
We have produced the same spin squeezed states via op- 
timal control, with small but significant reductions in 
both squeezing and yield. This suggests that gains from 
reduced decoherence (optimal control is as much as five 
times faster) is offset by increased sensitivity to control 
errors. 

In conclusion we have implemented a scheme for opti- 
mal control of the hyperfine spin-angular momentum of 
Cs atoms in the F = 3 ground state. Control Hamilto- 
nians were designed to produce a range of target states, 
applied in the laboratory and evaluated by measuring 
the resulting density matrices. Typical yields fall in the 
range 0.8-0.9. Among the states produced were a number 
of spin squeezed states, which allowed direct comparison 
of our optimal control approach to an adiabatic scheme 
that is more robust to errors in the control fields. In fu- 
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FIG. 4: (Color online) Spin squeezing by adiabatic control, 
(a) Normalized squeezing parameter vs. final magnetic field 
for the squeezed and anti-squeezed components. Dashed lines: 
perfect squeezing. Open symbols: predictions from our the- 
oretical model. Filled symbols: experimental results, (b) 
Target and (c) measured Wigner functions corresponding to 
the smallest observed £. 



ture experiments we plan to use a combination of rf and 
microwave magnetic fields to control the entire 16 di- 
mensional state space for the Cs 6Si/2(F = 3,4) ground 
manifold. Preliminary studies indicate that this system 
is fully controllable on a timescale of a few tens of mi- 
croseconds with easily available control fields. This will 
provide an important tool for the encoding and manipu- 
lation of qubits and qudits embedded in a larger atomic 
ground manifold. In the longer term it is also interesting 
to consider if control Hamiltonians of the form used here 
can be achieved for collective spins, for example through 
coherent optical feedback [2(| or through atom-atom in- 
teractions in a quantum-degenerate gas [21 1 . 
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